Magnetism without magnetic impurities in oxides ZrQ 2 and Ti02 
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We perform a theoretical study of the magnetism induced in transition metal dioxides ZrC>2 and 
Ti02 by substitution of the cation by a vacancy or an impurity from the groups 1A or 2A of the 
periodic table, where the impurity is either K or Ca. In the present study both supercell and 
embedded cluster methods are used. It is demonstrated that the vacancy and the K-impurity leads 
to a robust induced magnetic moment on the surrounding O-atoms for both the cubic Zr02 and 
rutile Ti02 host crystals. On the other hand it is shown that Ca-impurity leads to a non magnetic 
state. The native O-vacancy does not induce a magnetic moment in the host dioxide crystal. 

PACS numbers: 77.80.Bh 75.30.Et 
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I. INTRODUCTION 

We are all used when we refer to magnetic materials 
to have in mind systems which at least contain the es- 
sential brick, the magnetic impurities such as Fe, Ni or 
Mn. However, it has been reported several times that 
magnetism can also be observed in materials which does 
not contain such magnetic impurities. Ferromagnetism 
in this new class of materials is referred to as or 
intrinsic ferromagnetism. Experimentally, it has been 
shown that thin films of HTO2 and ZrO^ 1 ^ exhibit very 
high Curie temperature (above room temperature). Un- 
expected ferromagnetism has been also reported in ir- 
radiated graphite 3 -^ and even in La-doped hexaborides 
such as CaB6 5 . Theoretically, ferromagnetism was pre- 
dicted in systems, where the Fermi level lies in a flat 
band§i Recently it has been shown that CaAs is a good 
candidate for flat band ferromagnetism^. 

However, the intrinsic ferromagnetism is often a sub- 
ject of controversy. For example, it has been argued that 
the ferromagnetism in CaBg is due to the presence of 
magnetic impurities such as Fe^. In addition, intrin- 
sic ferromagnetism often appear difficult to be stabilized 
and controlled as it strongly depends on the sample his- 
tory: materials grown under similar conditions may lead 
to non-magnetic samples. Until now, there is no clear 
understanding of the underlying mechanism which is at 
the origin of the observed magnetism. From a theoret- 
ical point of view, it has been shown that point defects 
as cation vacancies may be at the origin of the mag- 
netism in these material a 10 : 11 ' 12 : 13 . The vacancy induces 
in these materials local magnetic moments on the neigh- 
bouring oxygen atoms which then interact with relatively 
extended exchange couplings. The oxygen vacancies were 
originally also suggested as a source of the magnetism in 



HfO^ 1 -, but the density functional calculations-^ found 
no evidence for the formation of magnetic moments 
around such defects. 

The long-range ferromagnetic order is possible only 
if the exchange couplings are extended enough (beyond 
the percolation threshold). Osorio-Guillen et al*£ have 
shown that in CaO oxide the ferromagnetism induced by 
cation vacancies is possible only if the density of vacan- 
cies is sufficiently large (above 5%). They found that the 
formation energy of a vacancy is relatively high, which 
implies that the cation vacancy induced ferromagnetism 
is unlikely. 

A progress has been done recently in the understand- 
ing of the ferromagnetism^ in the framework of a 
minimal model which includes on equal footing the effects 
of disorder and electron-electron correlation effects. The 
effect of the cation vacancy or non-magnetic impurity on 
neighbouring oxygen orbitals (the correlated disorder) is 
the relevant part of the theory. This study has predicted 
that very high Curie temperatures can be reached for 
realistic parameters. It was shown that the optimal situ- 
ation for ferromagnetism depends mainly on two param- 
eters (i) the position of the impurity band which should 
be on the edge to split from the valence band for the 
hole doped case, and (ii) the density of holes per defect 
should be of the order of three. In order to circumvent 
the difficulties of controlling intrinsic defects such as va- 
cancies, an alternative way which consists in the substitu- 
tion of the four-valent cation A 4+ in dioxides such as AO2 
(A=Ti, Zr, Hf) by a single-valent cation of the group 1A 
of the periodic table has been proposed. The substitution 
of the four-valent host A A+ by a single-valent impurity 
X + , where X= K, Li, Na, or Cs, introduces three holes 
per impurity. However, we should emphasize that it is 
not granted that the system will become magnetic and 



2 



that the impurity band will be located near the top of 
the valence band. 

In this paper we shall address a question concerning 
the formation of magnetic moments in the host dioxides 
doped by non-magnetic impurities. To this end we per- 
form an ab initio study of the effect of substitution of 
the cation by an element of the group 1A or 2A. We 
shall employ two different approaches, both based on the 
local density approximation (LDA) to the density func- 
tional approach: (i) the supercell approach using the 
full-potential linearised augmented plane-wave method 
(FLAPW)^£, and (ii) the embedded cluster method 
(ECM}H based on the Green function approach as im- 
plemented recently in the framework of the tight-binding 
linear muffin-tin orbital (TB-LMTO) method 21 . The sin- 
gle impurity in the supercell with increasing sizes corre- 
sponds to decreasing but still finite impurity concentra- 
tions. On the other hand, the embedded cluster approach 
allows us to address the problem of moment formation 
for a single impurity in the infinite host crystal. Both 
approaches, however, allow for electronic relaxations on 
the sites neighboring impurity which is an essential re- 
quirement for the formation of induced non-local mag- 
netic moments such as the A-vacancy in AO and AO2 
oxide o n ' 12 i 13 . We have also investigated the case of two 
nearest-neighbor impurities that have a common nearest 
neighbor oxygen site and the effect of group 2A impu- 
rities (such as Ca) which introduce only two holes per 
Ca-impurity when substituting the A-cation in AO2 ox- 
ides. 



II. FORMALISM 

Below we briefly summarize the supercell and ECM 
approaches used in the present paper. 

A. The supercell approach 

We consider a tetragonal supercell consisting of two 
conventional cells for the fluorite structure of Zr02 and 
four tetragonal units for the Ti02 rutile structure. It is 
constructed by doubling the area of the basis in the xy 
plane, and in the case of Ti02 also along the c axis. The 
supercell contains 8 molecular units of the host dioxide 
and its lattice parameter are (y/2a, V2a, a) for the fluo- 
rite structure and (v2a, y/2a, 2c) for the rutile structure. 
The substitutional impurity is introduced in the centre of 
the supercell XA7O16 which corresponds to an effective 
impurity concentration of 12.5%. We neglect the change 
of the lattice constant due to impurity as well as local 
relaxations. We have used the experimental lattice con- 
stants a=5.256 A for fluorite ZrO^ and a = b = 4.59 A, 
c = 2.96 A for the rutile TiO^. We have also con- 
sidered larger supercells corresponding to lower effective 
concentrations, namely, the XA31O64 supercell (concen- 
tration of 3.125%) with lattice parameters (2a, 2a, 2a) 



for the fluorite and (2a, 2a, 4c) for the rutile structures, 
respectively. 

The spin-polarized electronic structure was obtained 
by the FLAPW method as implemented in the WIEN2k 
program package^ with the Perdew-Wang exchange- 
correlation potential-^. We employed a plane- wave cutoff 
energy E cut = 16 Ry and 40 special points in the irre- 
ducible wedge of the Brillouin zone of the large super- 
cell with 96 atoms and correspondingly higher number of 
points for smaller supercells. 

The local magnetic moments are obtained by inter- 
grating over the muffin-tin spheres. Radii 2.35 a.u. and 
1.7 a.u. for K/Zr and O in Zr02, respectively and 2.0 a.u. 
and 1.6 a.u. for K/Ti and O in Ti02, respectively were 
used. 



B. The embedded cluster approach 

The embedded cluster method (ECM) is based on the 
Green function approach. The solid is divided into a clus- 
ter containing several tens or hundreds of atoms and the 
rest of the system. The cluster consists of the central (im- 
purity) atom (or central atoms, e.g., pair of impurities) 
and their neighbors whose potentials are allowed to differ 
from the host atoms due to the presence of perturbation. 
It is useful to make calculations for a sequence of grad- 
ually growing clusters in order to check how calculated 
quantities converge to a well-defined value. 

We first find the self-consistent solution for an ideal 
translationally invariant host as described by the LDA- 
Hamiltonian and corresponding Green function 

G^°\ The perturbation Vpp is localized, limited to the 
cluster region characterized by the projection operator 
P on the states of the cluster. For the selfconsistent 
LDA-solution of the perturbed problem we need only the 
PP-block of its Green function G described by the Hamil- 
tonian H = H° + Vpp, namely, 

Gp P = [(G^r 1 -Vp P ] PP (1) 

which is evaluated in real space from the known values of 
Gp P calculated directly in the reciprocal space. As the 
atomic charges in the cluster region are different from 
those in the host, we have to include the necessary cor- 
rections to the host Madelung potentials. We have im- 
plemented the embedded cluster method within the TB- 
LMTO scheme 2 ^. 



III. RESULTS AND DISCUSSIONS 

We have considered two oxides, Zr02 in its cubic (flu- 
orite) phase and Ti02 in the rutile structure. Two kinds 
of impurities were considered, namely vacancies and K- 
dopants in both cases. Test calculations were also done 
for Ca-dopants in Ti02, and for the nearest neighbor 
X— O— X complex in Zr02 (X=vacancy, K-impurity). 
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The results are summarized in Tables |TT] and IIVI for 
Zr0 2 and for Ti0 2 , respectively, using the FLAPW 
method for two effective impurity concentrations of 
3.125% and 12.5% corresponding to 96- and 24-atom su- 
percells. In Table Hill we present the results for a single 
impurity in the cubic Zr0 2 as calculated by the ECM 
approach with the charge selfconsistency for 59 sites of 
the embedded cluster (see Table [I]). 

TABLE I: Coordination spheres around the impurity X re- 
placing Zr atom located at the origin in Zr02. The symbols 
E denote empty spheres used in the TB-LMTO calculations. 
We show the coordinates r of only one of the atoms in the co- 
ordination sphere as the other follow from the symmetry 
of the cluster (multiplicity n). 



site 


X 





E 


Zr 


O 


E 


n 


1 


8 


6 


12 


24 


8 


r 


(000) 


|(111) 


|(200) 


f (220) 


f(311) 


f (222) 



In both approaches we have found no induced magnetic 
moment in the host due to native O-vacancies either in 
Zr0 2 or in Ti0 2 . This is in agreement with results of 
Das Pemmaraju and Sanvitc-^ for O-vacancies in Hf0 2 . 



A. Impurities in ZrC>2 

As seen from Table HH both vacancy and K-impurity 
on Zr-cation sublattice give rise to a well defined robust 
magnetic moment induced on neighboring shells of O- 
atoms. Note that the value of this moment could vary 
with the impurity concentration. We first discuss the case 
of a vacancy. For relatively low concentration of 3.125% 
the total induced moment in the supercell is 4 [Ib- We 
found that the induced moment decreases only slightly 
to 3.86 [Ib for a higher concentration of 12.5%. The rel- 
evant part of the induced magnetic moment is mainly 
located on the first shell of host O-atoms. The induced 
moments on the next oxygen shell are an order of mag- 
nitude smaller. The moment on Zr-atoms neighbouring 
the defect are negligible. Despite the smallness of induced 
moments on O-sites, the total moment is large because of 
the number of equivalent atoms in each shell neighboring 
the impurity is large (8 in the first shell, 24 in the second 
shell). These results are in a good agreement with those 
of Ref. UU for vacancies in the monoclinic Hf0 2 crystal. 
Of course, the magnitude of the induced moments could 
be quantitatively influenced by possible lattice relaxation 
as discussed in Ref. 

We now discuss the case of K-cation impurity. In this 
case we have also found an induced moment but slightly 
smaller than that of vacancy. Here, on the contrary, the 
system tends to a halfmetallic behavior for sufficiently 
high concentration of defects. The total induced mo- 
ment slightly increases with the impurity concentration. 
As expected the total induced moment is smaller roughly 
by one hb smaller as compared to that of vacancy. In- 
deed, K-impurity brings only three holes as compared 



TABLE II: The total magnetic moment per cell (M ce ii) and 
local magnetic moments on the impurity (mx) and on the 
nearest-neighbor host atoms adjoining it for two cubic su- 
percells, XZrrOie and XZr3i064, corresponding to effective 
impurity concentrations of 12.5% and 3.125%, respectively. 
Here moi, mo2, and mz r i denote atoms in corresponding 
nearest-neighbor shells of the host. The impurity atoms are 
X= vacancy and K. All magnetic moments are given in units 

Of flB- 



impurity X 


m x 


moi 


mo 2 


mzri 


M C ell 


vacancy ( 12.5%) 




0.40 


0.05 


-0.03 


3.81 


vacancy (3.125%) 




0.38 


0.03 


-0.01 


4.00 


K ( 12.5%) 


0.10 


0.30 


0.05 


-0.02 


2.99 


K (3.125%) 


0.08 


0.21 


0.03 


-0.01 


2.67 



to four holes in the case of a vacancy. Similarly to the 
case of vacancy, magnetic moments induced on more dis- 
tant O-atoms are an order of magnitude smaller but they 
are nevertheless important as they help to overcome the 
magnetic percolation threshold. The magnetic moments 
on other sites (impurity and neighboring Zr atoms) are 
again negligible. Because the induced moment around 
the cation impurities is rather extended, the correspond- 
ing percolation threshold is expected to be much smaller 
as compared to that of conventional local magnetic mo- 
ments e.g. in GaAs:Mn or GaN:Mn alloys. This is a par- 
ticularly interesting result for the substitutional K-cation 
impurity which can be more easily controlled experimen- 
tally than the intrinsic defects. 

It is interesting to compare the results of supercell cal- 
culations corresponding to small but finite impurity con- 
centrations with those obtained for a single impurity by 
the ECM approach (see Table [TTT|) . The total induced 
magnetic moment in the insulator should be, strictly 
speaking, 4 /is and 3/is for vacancy and K-impurity, re- 
spectively. This is fulfilled very well in the case of vacancy 
but not so accurately for K-impurity. This indicates that 
a small part of the induced moment resides outside the 
chosen cluster. We note that the results of the ECM cal- 
culations could be also interpreted as the condition for 
the formation of a local induced moment in the system 
(the single- impurity limit). 

In the framework of the supercell approach we have 
also investigated the possibility of induced magnetic mo- 
ments due to Ca-substitutional impurity in Zr0 2 . We 
have found no induced moment in this case. This re- 
sult is in agreement with a general conclusion made in 15 , 
namely that the optimal situation for magnetism in diox- 
ides requires impurities which release three holes per de- 
fect into the valence/impurity band. Thus, we find a de- 
creasing tendency to form an induced magnetic moment 
in Zr0 2 if Zr is substituted by vacancy, K, and Ca. This 
result can be understood qualitatively as a consequence 
of decreasing the charge difference between the host and 
impurity cation. Such a qualitative explanation could be 
supported by a simple Stoner-picture model. We present 
in Fig. Q] the total densities of states (DOS) of the non- 
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FIG. 1: The total densities of states per formula unit for 
various impurities in the cubic nonmagnetic Zr02: (a) ZrC>2 
host; (b) vacancy in ZrCh; (c) K in Zr02; and (d) Ca in Zr02- 
Calculations were performed for the supercell XZris032 cor- 
responding to the effective impurity concentration of 6.25%. 
The energy zero indicates the Fermi energy. 



magnetic phase calcualted using the supercell approach 
in the limit corresponding to the effective impurity con- 
centration of 6.25% (XZri 5 032, X=vacancy, K, Ca). A 
strong reduction of the total density of states at the Fermi 
energy for Ca-impurity as compared to K-impurity and 
vacancy is clearly seen. This explains the decreasing ten- 
dency to form an induced magnetic moment. 

Using the ECM approach we have also investigated the 
case of two nearest-neighbor vacancies/K-impurities in 
the ZrC>2 host that have two equivalent nearest-neighbor 
oxygen sites. The total induced moment on O-sites in- 
side the cluster is strongly enhanced when compared to 
the case of a single impurity: 7.1 fiB for two vacancies 
and 4.85 \ib for two K atoms. This should be compared 
with the values of 3.93 hb and 2.62 hb found for a single 
vacancy and K-impurity (see Table IIII[) . The local mo- 
ment on O atoms neighboring the impurity is strongly 
enhanced, namely 1.22 (is and 0.77 hb for the vacancy 
and K-impurity, respectively. The other induced mo- 
ments are similar to those found for a single impurity 
case. One could conclude that the nearest-neighbor X- 
O-X cluster (X=vacancy or K) induces in Zr02 host a 
larger and more extended magnetic moment. 



TABLE III: The total magnetic moment inside the cluster 
M c i ust and local magnetic moments on the impurity (mx) 
and on the nearest neighbor host atoms adjoining it for the 
cubic XZr02 crystal with a single impurity calculated using 
the embedded cluster method. Here moi, mo2, and mzri 
denote atoms in corresponding nearest-neighbor shells of the 
host. The impurity atoms are X= vacancy, and K. All mag- 
netic moments are given in units of /j,b- 



impurity X 


m x 


moi 


mo 2 


m Zri 


Mclust 


vacancy 


0.005 


0.421 


0.027 


-0.007 


3.927 


K 


0.158 


0.244 


0.023 


-0.005 


2.616 
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FIG. 2: (a) The total density of states per formula unit for 
Zr02 host and (b) the spin resolved densities of states fro 
KZrrOie supercell. (c) The spin-resolved local densities of 
states per formula unit on the first O-shell surrounding the 
impurity for the cubic KZr70i6 supercell. The effective K- 
impurity concentration is 12.5%. The energy zero indicates 
the Fermi energy. 



The spin resolved total DOS and local DOS on O- 
sites for substitutional K-impurity in Zr02 are shown in 
Fig. for the KZryOi6 supercell (12.5% of impurities). 
The substitutional K-impurities in Zr02 polarize the host 
band giving rise to well split spin-subbands. The induced 
impurity peak close to the valence band edge is seen in 
the minority DOS while the majority DOS is influenced 
by impurities only negligibly (compare Fig. [5^, andf^b). 
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This effect is even more clearly seen in the O-local DOS 
shown in Fig. [2b from which it is found that the valence 
band is formed predominatly by the oxygen p-states. The 
above results are again in a qualitative agreement with 
those of Ref. [l5|. Namely, the presence of an impurity 
band close to the top of the valence band is a precursor 
for the appearence of magnetism in dioxides. This is the 
case for both a vacancy and substitutional K-impurity 
(see Fig. [Tb^c). In contrast, for Ca substitution, the im- 
purity band is shifted deeper into the valence band, which 
prevents the formation of a magnetic moment. 



B. Impurities in TiC>2 

We have also investigated the possibility for non- 
magnetic impurity induced magnetism in TiC>2 . We used 
similar supercells, namely XTiyOi6 an d XTi3iC>64 as in 
the case of ZrC>2 (effective concentrations are 12.5% and 
3.125%). The main results (see Table [TV} are in qual- 
itative agreement with those for the cubic ZrC>2 struc- 
ture. Nevertheless, for higher concentration of vacancies 
(12.5%) no magnetic state has been found in TiC>2. For 
smaller concentrations both vacancy and substitutional 
K-impurity induce magnetic moments on the surround- 
ing O-sites and the size of total induced moments are 
similar in ZrC>2 and TiC>2. There are, however, some 
quantitative differences due to different lattices. In the 
rutile structure there are four equivalent nearest-neighbor 
O-sites and two other next nearest-neighbor O-sites with 
almost the same distance from the Ti- /impurity site. 
Lattice anisotropy leads to a different distribution of the 
induced moments. 

TABLE IV: The total magnetic moment per cell (M ce ii) and 
local magnetic moments on the impurity (mx) and on the 
nearest neighbors host atoms adjoining it for two cubic su- 
percells, XTLOie and XTi3iOe4, corresponding to effective 
impurity concentrations of 12.5% and 3.125%, respectively. 
Here moi, mo2, and mTi ± denote atoms in corresponding 
nearest-neighbor shells of the host. The impurity atoms are 
X=vacancy, and K. All magnetic moments are given in units 
of hb- 



impurity X 


m x 


moj 


mo 2 


niTij 


M ce ll 


vacancy ( 12.5%) 




0.00 


0.00 


0.00 


0.00 


vacancy (3.125%) 




0.44 


0.28 


-0.04 


3.46 


K ( 12.5%) 


0.27 


0.05 


0.32 


0.00 


2.10 


K (3.125%) 


0.35 


0.31 


0.18 


0.01 


2.91 



We have also checked the effect of local structure relax- 
ations in the vicinity of impurity for the case of KTi70i6- 
The distance between the K-impurity and four near- 
est neighbor 01 oxygen sites increases from 1.95 A to 
2.27 A and between two 02 oxygen sites from 1.98 A to 
2.15 A . The existence of induced moment due to K- 
impurity is not influenced, the total induced moment 
changes from 2.10 [is to 2.08 fiB, the local moments on 
Ol-sites change from 0.32 {1b to 0.26 [is while the local 



moments on 02-sites remain unchanged (0.05 (Us). 

The induced local moments are different for different 
impurity concentrations, e.g. on 02 sites is local mag- 
netic moment 0.19 [is for a concentration of 3.125% but 
becomes —0.05 \ib when the impurity concentration is 
12.5%. The reason for this variation is not clear to 
us. For example, one could speculate that this could 
be a result of different hole concentrations in both spin- 
subbands. We should also mention a relatively large local 
moment on K-impurity which appear to be comparable 
to that induced on individual O-sites. It is few times 
larger than the local K- moment in Zr02 (see Table Hi]) . 

The total and local DOS for KTiyOig supercell are 
shown in Fig. [3J The general trend is similar to that 
found for K-impurities in Zr02: there is an impurity 
band in the minority DOS but less pronounced thus indi- 
cating a slightly weaker disorder due to K-impurities in 
Ti02 as compared to Zr02 (in terms of model language 
of Ref.QJ). The majority band is again only weakly per- 
turbed and the system is metallic with the Fermi energy 
lying inside the valence band for both the majority and 
minority bands. 

The spin resolved local DOS on 02 is weakly spin 
dependent for KTi70i6 supercell (see Fig. [5b). This 
results in a very small induced moment on these sites 
(—0.05 ixb)- On the other hand, the 01 local DOS shows 
a halfmetallic behavior and gives rise to a larger induced 
moment (0.28 /is). For lower concentration we observe a 
halfmetallic behavior for both 01 and 02 spin-resolved 
DOS (see Fig. [3Ji). This results in a similar size for the 
moment on 01 atoms but a significantly increased 02 
moment (see Table HVf . 



C. Local magnetization around impurities in ZrC>2 

The local magnetization m(r) in [110]-plane around 
the impurity in real space is shown in Fig.[4]for KZryOi6 
supercell. The small local moment on K is in the centre 
of the plot with strongly damped Friedel-like oscillations. 
Such a damping is due to halfmetallic behavior of this 
system, i.e., the position of the Fermi energy lies in the 
gap of the majority band. This feature is similar to 
that observed in the magnetic exchange couplings in the 
diluted magnetic semiconductors 22 . Large increase of 
the magnetization density in the neighborhood of the 
O-sites indicates the existence of an induced moment, 
few times larger than the local K-moment (see Table [TTJ) . 
The induced moments on the second shell of O-sites are 
also shown. Their magnitude is small but because the 
number of sites in the shell is relatively high (24), their 
contribution is non-negligible. 
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FIG. 3: The total densities of states per formula unit for (a) 
TiC>2 host with the rutile structure and (b) spin resolved den- 
sity of states for KTi70i6- The local spin-resolved densities of 
states on the first two O-shells surrounding the impurity in the 
ferromagnetic cubic (c) KTiyOie and (d) KTiis036 that cor- 
respond to concentrations of 12.5% and 6.25% Respectively. 
The energy zero indicates the Fermi energy. 

IV. CONCLUSIONS 

In conclusion, we have shown successfully that a new 
path to controlled d°-fcrromagnetism can be achieved 
by cation doping of oxides with well defined elements of 
the periodic table. Specifically, we have shown that the 
promissing systems are both cubic ZrC>2 and rutile T1O2 



materials doped with K (group 1A). The methods which 
have been used are both the supercell approach and the 
embedded cluster method. The main conclusions of our 
study are: (i) The magnetic moment is mainly induced on 
the first oxygen shell neighbouring the defect, (ii) A reso- 
nance peak near the top of the valence band in the minor- 
ity channel is found, while the majority channel remains 




FIG. 4: The local magnetization for the ferromagnetic cubic 
KZrrOie supercell in the [110]-plane. The K-impurity is in 
the center and two nearest neighbor O-atoms in the first O- 
shell, while four O-atoms in the second O-shell are visible at 
[Oil] edges. 



mainly unperturbed. This result is in agreement with a 
recent theoretical study based on model calculations 15 , 
(iii) The substitution by Ca (group 2A) has lead to a 
non magnetic ground state, (iv) We have verified for 
a specific case that the total induced moment is robust 
with respect to lattice relaxations, (v) Our results also 
indicate that these effects are not inherent to a specific 
lattice structure. 

The interesting and crucial issues which have to be ad- 
dressed in the future are the impurity formation energies 
and the exchange integrals between induced moments. 

This work has been done within the project AVOZ1- 
010-0520 of the AS CR. We acknowledge the support 
from the Grant Agency of the Czech Republic, Contract 
No. 202/07/0456, COST P19-OC150 project, and sup- 
port from the Grant Agency of the Academy Sciences of 
the Czech Republic (A100100616). 
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